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■ New local "hybrid" functionals proposed by V. V. Karasiev in [J. Chem. Phys. 118, 8567 (2003)] 
I are benchmarked against nonlocal hybrid functionals. Their performance is tested on the total and 

'/^ ■ high occupied orbital energies, as well as the electric moments of selected diatomic molecules. The 

(N- new functionals, along with the Hartree-Fock and non-hybrid functionals, are employed for finite- 

' difference calculations, which are basis-independent. Basis set errors in the total energy and electric 

^ , moments are calculated for the 6-311G, 6-311G-f -f G(3df,3pd) and AUG-cc-pVnZ (n=3,4,6) basis 

^ I sets used in conjunction with the Hartree-Fock and conventional density functional methods. A 

I comparison between the results of the finite-difference local "hybrid" and basis set nonlocal hybrid 

S , functional shows that total energies of local and nonlocal hybrid functionals agree to within the 

^ ■ basis set error. Discrepancies for multipole moments are larger in magnitude when compared to the 

(-H ■ basis set errors, but still reasonably small (smaller than errors produced by the 6-311G basis set). 

^ , Thus, we recommend using the new local "hybrid" functionals whenever the accuracy is expected 
to be sufficient, because they require a solution of just differential Kohn-Sham equations, instead of 

^ ■ integro-differential ones in the case of hybrid functionals. 
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I. INTRODUCTION 

> 

_ The investigation of electronic structure by means of theory represents an important and highly active area of 
QQ ■ research. Calculations of electronic structures are frequently performed by means of the density functional theory 
\^ ' (DFT)ji"— DFT in the form of the Kohn-Sham (KS) method provides an approach to treat the Schrodinger equation 
within an exact formalism that bypasses numerically prohibitive calculations based on many-electron wavefunctions. 
Within DFT, exchange and correlation effects are accounted for by means of energy functionals of the electron 
ps) i density - the exchange-correlation functionals, and their functional derivatives with respect to the electron density - 
T-H ■ the exchange-correlation potentials. These functionals are approximated in standard DFT methods. 
^ ] Hybrid functionals which combine the nonlocal Hartree-Fock exchange with local KS exchange-correlation potentials 
• ^ ■ provide a significant improvement over conventional DFT approaches with respect to binding energies and other 
■ properties. Unfortunately, such hybrid methods no longer fit the framework of the DFT formalism due to presence 
of a non-multiplicative operator in corersponding one-electron equations. Moreover, the incorporation of a nonlocal 
Hartee-Fock exchange operator, defined by its action on a particular electron orbital as 



/■ ^*(r')V'»(r') , 

-^5(a.,a,)^,(r) / ^J^dv' , (1) 
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converts the system of differential KS-DFT equations into a system of integro- differential ones. The solution of such a 
system of integro-differential equations is more tedious procedure compared to a solution of KS-DFT equations with 
local multiplicative potential, which is the same for all KS orbitals. E^^ in Eq. ([T]) is the exact exchange energy (the 
exchange energy of KS determinant) defined by the Hartree-Fock expression: 

1 " 

; — ardr . (2) 

|r — r 

To alleviate the non-locality of the HF exchange operator, the exact-expression approximate exchange (EFAX) 
methods were developed. They all share the following common properties)^ (i) the exchange energy expression is the 
exact one defined by Eq. ([5]), and (ii) corresponding exchange potential is a local, multiplicative operator. Among the 
first EFAX approximations, we should mention method introduced by Krieger, Li, and lafrate (KLI),^ the localized 
Hartree-Fock (LHF) method and the common energy denominator approach (CEDA) introduced by Gorlingl and 
Gritsenko and Baerends^ respectively. Later, self-consistent a (SCa)^"— and asymptotically-adjusted self-consistent 
a (AASCa)^ local exchange functionals were developed. Since the exchange energy in all of these methods is the 
exact one, defined by Eq. the self-interaction error is eliminated, while multiplicative (or local in that sense) 
exchange potential guarantees that these methods are fully in the framework of DFT formalism. Another important 
aspect is that the EFAX methods satisfy the Levy-Perdew virial relationj ^'^'^^ 

A solution to the problem of non-local operator in hybrid functionals was proposed and implemented for the first 
time in Refi^, when new local "hybrid" exchange-correlation functionals obtained by replacement of the nonlocal 
Hartree-Fock (HF) exchange in conventional hybrid models by one of the EFAX methods (the SCa and the AASCa 
approximations are used in Refi^). Such functionals are fully within the framework of DFT formalism because 
the functional derivative of EFAX functionals w.r.t. electronic density is a local multiplicative potential. From 
the computational point of view, the local "hybrid" functionals can be viewed as more efficient alternative to the 
conventional hybrid methods. 

The same scheme was later used to constract local "hybrid" functionals based on replacement the HF exchange in 
non-local hybrids by the LHF approximation with local exchange potential. ^^^^^ 

A comparison between non-local and local hybrids based on the SCa and AASCa methods which has involved 
total and HOMO energies, vibrational frequencies and bond lengths was made for diatomic molecules in Ref.— . The 
comparison of other properties, in particular electric moments can afford a more detailed measure of the closeness of 
local and nonlocal hybrid functional performance. 

In the present work we apply the local "hybrid" functionals introduced in^ for finite-difference (FD) calculations of 
total and HOMO energies, and dipole and quadrupole electric moments for selected diatomic molecules. To compare 
the local "hybrid" functional results with corresponding nonlocal hybrid calculations, for which the basis set (BS) 
technique is used, an estimate of basis set error (BSE) for each magnitude compared is needed. For this purpose BSEs 
are calculated as the difference between the FD and BS results for the HF and conventional DFT exchange-correlation 
functionals. 

In Sec. II the local "hybrid" functionals used in calculations are introduced and definitions of the calculated 
quantities are established. Numerical results are discussed in Sec. III. In Sec. HI A the BSEs for total energy 
and electric moments of standard basis sets available in GAUSSIAN packag e^^i^^ are presented. Local "hybrid" 
functional self-consistent results for total energy, multipole moments and one-electron energies are compared to the 
values obtained from calculations with nonlocal ones in Sec. HI B. Finally, our conclusions are presented in Sec. IV. 

II. METHOD 

Two nonlocal hybrid functionals, PBFO^S based on the Perdew-Burke-Ernzerhof (PBE) DFT exchange-correlation 
functional, and one parameter Becke-Lee-Yang and Parr hybrid, BILYP)^"— are used to construct local "hybrids". 
The nonlocal HF exchange in these non-local hybrids was replaced by one of the EFAX term. First, by SCa-PBF 
and SCa-B88 functionals respectively (see Refj^i for details), and second, by AASCa exchange (AA-m2 model from 
Ref.-'^^-.). For example, local SC-BILYP "hybrid" obtained from non-local BILYP functional by replacement of the 
non-local exchange term by the local SCa-B88 functional has the following form 




iSC-BlLYP _ 

'xc — 



(3) 
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with ao = 0.25. Expression of E^'~^°' exchange functional is exact one defined by Eq. ^ , while exchange-correlation 
potential corresponding to functional of Eq. is a local multiplicative operator 



V. 



LSD I A„,LYP 



A.-^ (4) 



where potential wSCa-B88 jg defined in Ref.-ii. On the whole, four local "hybrids", SC-PBEO, AAni2-PBE0 and SC- 
BILYP, AAm2-BlLYP were employed for numerical calculations of total and HOMO energies and electric moments 
of selected diatomic molecules. 

The multipole moments were calculated according to the following definition (see Refsi^i^). Equation for the 
independent component for each electric moment for a system with axial symmetry 

MC^) = ^A^A. - f r^Pk{z/r)p{v)dPv, (5) 

A 

where Pk are Legendre polynomials of degree k, was used in finite-difference code to calculate the first two moments 

M^^^=li^fx,, (6) 

and 

= e = e,,. (7) 

The multipole moments were determined from the basis set calculations, 

fia^Y ZaRac ~ / rapir)(fr, (8) 



and 



Qafs = Y^ZARAaRAp - j rari3p{T)Sr, {a, (3 = X , y , z) , (9) 



A 

They can be related to components of Eqs. ©-j!]) evaluated by the FD code (see Ref.--^): 

= Qzz - ^(Qxx + Qyy). (10) 

Relations for higher moments can be found in Ref.^^. In Eqs. ([5]), ([H])-®, Za is the charge of nucleus A and ii^'s are 
its Cartesian coordinates and r^'s are {x,y,z) coordinates. 

The moments Af*^^^ and Af^^^ in finite-difference calculations and fiz, ^zz in basis set calculations were determined 
with respect to the molecular center of mass (GAUSSIAN keyword Symmetry=CenterOfMass). 



III. RESULTS AND DISCUSSION 



A. Basis set error for the HF and DFT methods 



Local "hybrid" exchange-correlation functionals discussed in the present work have been implemented in the spin- 
restricted FD full-numerical program by Kobus, Laaksonen and Sundholm^ while the non-local hybrid functional 
calculations were performed using basis set based GAUSSIAN packageJ^ii^ The accuracy achieved in FD calculations 
is about few /iHartree. The error of basis set truncation is usually higher. 

On Figure [T] the total energies of finite-difference calculations for the HF method and for the PBE and BLYP 
functional DFT calculations are presented. These values are compared to the basis set results for different basis sets. 
Corresponding BSEs defined as difference between FD and BS total energies are plotted. As one would expect, these 
values, due to the variational principle, are always positive (a numerical error of the FD code is smaller than the basis 
set truncation error). Values corresponding to the 6-311G basis set are quite similar for the HF and DFT methods 
and they are relatively small for the H2 and BH molecules (an order of 0.005-0.01 Hartrees). These errors increase for 
other systems up to 0.05-0.1 Hartrees. Augmentation of 6-311G basis with polarization and diffuse functions decreases 
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TABLE I: Orbital energies {—Si) of the CO molecule from the finite-difference HF and PBE, BLYP DFT calculations and 
BSE" for different basis set calculations (in Hartrees). 



system method 


ef*7BSE 


ef^^/BSE 


ef^^^/BSE 


la 


FD: [169xl93;40] 


20.6645 


18.8576 


18.8807 




6-311G 


-0.0192 


-0.0045 


-0.0034 




6-311++G(3df,3pd) 


-0.0032 


0.0031 


-f-0.0028 




AUG-cc-pV6Z 


0.0000 


0.0002 


-H0.0002 


2a 


FD: [169xl93;40] 


11.3600 


9.9997 


10.0249 




6-31 IG 


-0.0320 


-0.0073 


-0.0062 




6-311-h-hG(3df,3pd) 


-0.0009 


0.0043 


+0.0038 




AUG-cc-pV6Z 


0.0000 


0.0001 


+0.0002 


3(7 


FD: [169xl93;40] 


1.5216 


1.0769 


1.0700 




6-311G 


-0.0427 


-0.0416 


-0.0395 




6-311-H-G(3df,3pd) 


-0.0013 


-0.0010 


-0.0010 




AUG-cc-pV6Z 


-1 ■ 10"^ 


5 • lO-** 


+1 ■ 10~^ 


4(7 


FD: [169xl93;40] 


0.8045 


0.5192 


0.5144 




6-31 IG 


4-0.0044 


0.0095 


+0.0111 




6-311-H-HG(3df,3pd) 


-0.0007 


-0.0005 


-0.0005 




AUG-cc-pV6Z 


-1 • 10"^ 


1 ■ 10"^ 


+1 ■ 10"^ 


Itt 


FD: [169xl93;40] 


U.64U0 


0.4JOU 


A Aon^ 
0.4JU ( 




6-31 IG 


-0.0059 


-0.0066 


-0.0056 




6-311++G(3df,3pd) 


-0.0005 


-0.0002 


-0.0003 




AUG-cc-pV6Z 


0.0000 


5 • 10-^ 


-0.0001 


5(7 


FD: [169xl93;40] 


0.5549 


0.3322 


0.3311 




6-311G 


-0.0023 


0.0059 


+0.0071 




6-31H-+G(3df,3pd) 


-0.0004 


-0.0001 


-0.0002 




AUG-cc-pV6Z 


0.0000 


8- 10"'' 


+1 ■ 10"^ 



"BSE stands for basis set error for orbital energy estimated by difference £i[BS] — £i[FD]. 

an error by the factor of 2 to 10 for the molecules presented on Figure [T] (factor of 2 for the F2 molecule and factor 
of 10 for the CN", CO, N0+ and N2 systems). 

Finally, the BSEs of the largest AUG-cc-pV6Z basis sets are in fourth decimal place in most cases (with few 
exceptions, where the value is smaller). The AUG-cc-pV6Z basis set error in DFT calculations for the HCl molecule 
is larger (0.0012 Hartrees). 

As one can see from the Figure [1] the errors in HF method is smaller than the errors in DFT calculations when the 
AUG-cc-pV(T,Q,6)Z basis sets are employed. 

The errors introduced by incomplete basis set in one-electron orbital energies are presented in Table U for CO 
molecule as an example. The BSEs for the 6-311G basis are in third or in second decimal place. The addition of 
polarization and diffuse functions decrease these errors significantly. The AUG-cc-pV6Z basis set provides practically 
exact values to within four decimal places. 

Convergence of the dipole moment values with the basis set (and corresponding BSEs) for three molecules, CO, BE 
and NO"*" , are shown in Table [Hi The basis set error decreases very fast with increasing the basis set size. For the 
largest available basis set the moments are reproduced quantitatively within four decimal places. The same trends 
are observed for the quadrupole moment values, shown in Table IHII The convergence is fast for all the molecules 
presented in Table Hill The error for the largest basis set employed (AUG-cc-pV6Z) is in fourth decimal. 

In Table [TVl we present the dipole and quadrupole moments from FD and largest BS calculations corresponding to 
the HF and the DFT methods. The BSEs are calculated as /u[BS] - m[FD], and e[BS] - e[FD]. The AUG-cc-pV6Z 
basis set errors for dipole moment /i in most cases are in fourth decimal except for the DFT calculations of the HCl 
molecule and the CN" anion, where the errors are 0.0012 D and 0.0049 D respectively. 

The quadrupole moment basis set errors are also reproduced within four decimal places with a few exceptions. DFT 
calculations for the CN~ ion shows the largest error (0.0138-0.0150 DA). 

Here we should emphasize that the numerical errors of the finite-difference code for the multipole moments are in 
the order of 1 ■ 10~^ D for the dipole moment n and in the order of 1 • 10~^ DA for the quadrupole moment O. These 
values are much smaller than differences presented in Table IIVI which, consequently, can be considered as basis set 
errors corresponding to the multipole moments /i and O. 
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TABLE II: Dipole moments for the CO, BE and NO"'' diatomics from the finite-differences and basis set HF and DFT 
calculations and BSE" for various basis sets. 



system 


numerical method 


HF 


PBE 


BLYP 










BSE 


fi BSE 


/i BSE 


CO 


BS: 


6-311G 


-0.4974 


-0.2335 


0.1381 -0.0864 


0.1049 -0.0806 




BS: 


6-311++G(3df,3pd) 


-0.2708 


-0.0069 


0.2213 -0.0032 


0.1811 -0.0044 




BS: 


AUG-cc-pVTZ 


-0.2659 


-0.0020 


0.2238 -0.0007 


0.1845 -0.0010 




BS: 


AUG-cc-pVQZ 


-0.2642 


-0.0003 


0.2243 -0.0002 


0.1853 -0.0002 




BS: 


AUG-cc-pV6Z 


-0.2639 


0.0000 


0.2246 4-0.0001 


0.1857 +0.0002 




FD: 


[169xl93;40] 


-0.2639 




0.2245 


0.1855 


BF 


BS: 


6-311G 


0.4338 


-0.4410 


0.8692 -0.2028 


0.8379 -0.1952 




BS: 


6-311++G(3df,3pd) 


0.8740 


-0.0008 


1.0674 -0.0046 


1.0305 -0.0026 




BS: 


AUG-cc-pVTZ 


0.8764 -f0.0016 


1.0758 +0.0038 


1.0349 +0.0018 




BS: 


AUG-cc-pVQZ 


0.8744 


-0.0003 


1.0727 -f0.0007 


1.0335 +0.0004 




BS: 


AUG-cc-pV6Z 


0.8748 


0.0000 


1 0700 -L-O nnno 

l.U t ZZ -hU.UUUZ 


1 .Uooo -hU.UUUO 




FD: 


[169xl93;40] 


0.8748 




1.0720 


1.0331 


N0+ 


BS: 


6-311G 


-0.8157 


-0.1750 


-0.4066 -0.0927 


-0.4214 -0.0925 




BS: 


6-311++G(3df,3pd) 


-0.6439 


-0.0032 


-0.3152 -0.0013 


-0.3306 -0.0017 




BS: 


AUG-cc-pVTZ 


-0.6388 -hO.0019 


-0.3121 +0.0018 


-0.3269 0.0020 




BS: 


AUG-cc-pVQZ 


-0.6405 -h0.0002 


-0.3136 +0.0003 


-0.3286 0.0003 




BS: 


AUG-cc-pV6Z 


-0.6408 


-0.0001 


-0.3139 0.0000 


-0.3288 +0.0001 




FD: 


[169xl93;40] 


-0.6407 




-0.3139 


-0.3289 



"BSE stands for basis set error for the moments estimated by differences /i[BS] — /^[FD]. 

TABLE III: Quadrupole moments for the CO, BF and NO^ diatomics from the finite-differences and basis set HF and DFT 
calculations and BSE" for various basis sets. 



system immerical method HF PBE BLYP 









e 


BSE 


e 


BSE 


e 


BSE 


CO 


BS: 


6-311G 


-2.8917 


-0.8332 


-2.6523 


-0.6389 


-2.6637 


-0.6165 




BS: 


6-311++G(3df,3pd) 


-2.1252 


-0.0667 


-2.0769 


-0.0635 


-2.1071 


-0.0599 




BS: 


AUG-cc-pVTZ 


-2.0799 


-0.0214 


-2.0370 


-0.0236 


-2.0698 


-0.0226 




BS: 


AUG-cc-pVQZ 


-2.0590 


-0.0005 


-2.0181 


-0.0047 


-2.0517 


-0.0045 




BS: 


AUG-cc-pV6Z 


-2.0588 


-0.0003 


-2.0138 


-0.0004 


-2.0478 


-0.0006 




FD: 


[169xl93;40] 


-2.0585 




-2.0134 




-2.0472 




BF 


BS: 


6-311G 


-4.8487 


-0.6031 


-3.8964 


-0.5340 


-3.9176 


-0.5072 




BS: 


6-311++G(3df,3pd) 


-4.2418 


0.0038 


-3.3513 


0.0111 


-3.4025 


0.0079 




BS: 


AUG-cc-pVTZ 


-4.2422 


0.0034 


-3.3689 


-0.0065 


-3.4178 


-0.0074 




BS: 


AUG-cc-pVQZ 


-4.2417 


0.0039 


-3.3662 


-0.0038 


-3.4137 


-0.0033 




BS: 


AUG-cc-pV6Z 


-4.2449 


0.0000 


-3.3627 


-0.0003 


-3.4107 


-0.0004 




FD: 


[169xl93;40] 


-4.2456 




-3.3624 




-3.4104 




N0+ 


BS: 


6-311G 


-0.2997 


-0.9906 


-0.2349 


-0.8152 


-0.2499 


-0.7993 




BS: 


6-311++G(3df,3pd) 


0.6984 


0.0075 


0.5793 


-0.0010 


0.5509 


0.0015 




BS: 


AUG-cc-pVTZ 


0.6926 


0.0017 


0.5769 


-0.0034 


0.5471 


-0.0023 




BS: 


AUG-cc-pVQZ 


0.6946 


0.0037 


0.5819 


0.0016 


0.5513 


0.0019 




BS: 


AUG-cc-pV6Z 


0.6908 


-0.0001 


0.5800 


-0.0003 


0.5492 


-0.0002 




FD: 


[169xl93;40] 


0.6909 




0.5803 




0.5494 





"BSE stands for basis set error for the moments estimated by differences ©[BS] — 0[FD]. 

B. Local "hybrid" functional calculations: comparison to non-local hybrids 

Total energies, dipole and quadrupole moments for 10 diatomic molecules obtained from different exchange- 
correlation local "hybrid" functional methods are presented in Tables |Vl I VII and IVIIl The differences between values 
obtained from nonlocal hybrid functional calculations and values corresponding to the local "hybrid" functional are 
given in parenthesis. For example, differences presented in parenthesis of the second and third column of Table fVl are 
A = £;pBEo[BS] - Esc-PBEo[FD] and A = £'pbeo[BS] - EAAm2-PBEo[FD] correspondingly. 

The BSEs of basis sets employed for the nonlocal hybrid functional calculations are given for comparison. These 
values are estimated by difference between corresponding magnitudes of the BS and the FD DFT calculations and 
they are taken from Figure [T] and Tables |lTl IIIIl and IIVI 
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TABLE IV: Multipole moments for selected diatomic molecules from the basis set HF and PBE, BLYP DFT calculations and 
comparison to the finite-difference results. 



system numerical method 


HF 


PBE 


BLYP 








e 




e 




e 


H2 


BS: AUG-cc-pV6Z 


0.0 


0.6640 


0.0 


0.5786 


0.0 


0.5658 




FD: [I69xI93;28] 


0.0 


0.6639 


0.0 


0.5785 


0.0 


0.5658 




BSE" 


0.0 +0.0001 


0.0 +0.0001 


0.0 


0.0000 


BH 


BS: AUG-cc-pV6Z 


-1.7301 


-3.6021 


-1.5261 


-3.2986 


-1.5120 


-3.3423 




FD: [I69xf93;40] 


-1.7301 


-3.6017 


-1.5259 


-3.2986 


-1.5118 


-3.3413 




BSE 


+0.0000 


-0.0004 


-0.0002 


0.0000 


-0.0002 


-0.0010 


FH 


BS: AUG-cc-pV6Z 


-1.9219 


2.3300 


-1.7457 


2.2564 


-1.7511 


2.2503 




FD: [I93x235;50] 


-1.9218 


2.3299 


-1.7458 


2.2566 


-1.7510 


2.2506 




BSE 


-0.0001 


+0.0001 


+0.0001 


-0.0002 


-0.0001 


-0.0003 


HCl 


BS: AUG-cc-pV6Z 


-1.1797 


3.7679 


-1.0695 


3.5819 


-1.0474 


3.5473 




FD: [295x295;45] 


-1.1790 


3.7664 


-1.0683 


3.5803 


-1.0463 


3.5453 




BSE 


-0.0007 +0.0015 


-0.0012 +0.0016 


-0.0011 


+0.0020 


CN" 


BS: AUG-cc-pV6Z 


0.3928 


-4.4613 


0.7138 


-5.0930 


0.6767 


-5.1122 




FD: [I69xl93;40] 


0.3929 


-4.4608 


0.7177 


-5.0792 


0.6816 


-5.0972 




BSE 


-0.0001 


-0.0005 


-0.0039 


-0.0138 


-0.0049 


-0.0150 


NO+ 


BS: AUG-cc-pV6Z 


-0.6408 


0.6908 


-0.3139 


0.5800 


-0.3288 


0.5492 




FD: [I69xl93;40] 


-0.6407 


0.6909 


-0.3139 


0.5803 


-0.3289 


0.5494 




BSE 


-0.0001 


-0.0001 


0.0000 


-0.0003 +0.0001 


-0.0002 


N2 


BS: AUG-cc-pV6Z 


0.0 


-1.2507 


0.0 


-1.5359 


0.0 


-1.5660 




FD: [f69xl93;40] 


0.0 


-1.2509 


0.0 


-1.5364 


0.0 


-1.5663 




BSE 


0.0 +0.0002 


0.0 +0.0005 


0.0 +0.0003 


F2 


BS: AUG-cc-pV6Z 


0.0 


0.6758 


0.0 


0.8994 


0.0 


0.8559 




FD: [I93x235;50] 


0.0 


0.6753 


0.0 


0.8994 


0.0 


0.8558 




BSE 


0.0 +0.0005 


0.0 


0.0000 


0.0 +0.0001 



"BSE stands for basis set error for the moments estimated by differences /i[BS] — A*[FD] and 0[BS] — 0[FD]. 

Comparison of A's with corresponding values of BSE shows that total energies obtained from the AAm2-PBE0 
local "hybrid" functional calculations are coincide with energies coming from the nonlocal PBEO hybrid functional 
within the magnitude of BSE practically for all cases with one exception for the F2 molecule. For the SC-PBEO 
local functional the values of A are larger than corresponding values of BSE's but they do not exceed 0.0022 Hartree. 
The average absolute values of A's are 0.0012 and 0.0004 Hartree for SC-PBEO and AAm2-PBE0 local functional 
correspondingly and the average BSE for the PBE DFT functional has the value of 0.0003 Hartee. 

The situation is quite similar for the SC-BILYP and AAm2-BlLYP local functionals: the average absolute values of 
A's for the AAm2-BlLYP functional are very close to the average BSE (0.0005 and 0.0003 Hartree correspondingly). 
The maximum value of A for the SC-BILYP functional is 0.0021 Hartree for the F2 diatomic. 

A different trend is observed for the electric multipole moment calculations. Values of A which correpond to 
difference between non-local hybrids and local functional values are approximately one order of magnitude larger 
than corresponding values of BSE. For example, the average BSE value of BLYP functional for dipole moment is 
0.0010 D, the average absolute values of A for the SC-BILYP and AAm2-BlLYP functionals (i.e. average absolute 
differences between dipole moment values corresponding to the non-local BILYP hybrid and to the local "hybrid" 
functionals) are much larger (0.0753 D for the SC-BILYP and 0.0366 D for the AAm2-BlLYP). Dipole moment mean 
absolute error (MAE) of local "hybrid" functionals based on the asymptotically-adjusted potential (AAm2-PBE0 and 
AAm2-BlLYP) is about twice smaller as compared to the functionals based on the self-consistent a method (SC- 
PBEO and SC-BILYP). This could be explained by the fact that the AAm2 model for the exchange potential much 
better mimics the non-local behavior of the non-local Hartree-Fock exchange, that appears in conventional hybrids 
(see details irt^iliii^). 

In spite of fact that values of A for multipole moments are larger than corresponding BSE values, the relative 
discrepancies between local and non-local hybrid functional values are reasonably small (order of few percents for 
both dipole and quadrupole moments) , with exceptions for CN~, CO and N0+ dipole moment calculations when 
SC-PBEO and SC-BILYP local functionals are used (when the relative error is between 10-60 %). Comparing MAE 
from Tables IVll and IVlIl with BSE from Tables |ll] and HIT] respectively, we conclude that the finite-difference method 
combined with local "hybrid" functional produces dipole and quadrupole moments with better accuracy, than the 
corresponding non-local hybrid functional does in combination with the 6-31 IG basis set and worse than that with 
6-311++G**. 
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6-31 1G 

6-311++G(3df,3pd) ^-1.13 
AUG-cc-pVTZ 
AUG-cc-pVQZ 
AUG-oc-pVeZ 

FD 




1.14 a> 



-1.15 



-1.16 ^ 



-1.17 LU 



-25.15 



-25.25 B 

Q 
Ll. 

LJJ 



-100.1 



-100.2 t a) 



-460.4 5 



-460.8 LJJ 
-92.3 



-92.4 „ 
ID 

-92.5 



■92.6 
-92.7 
-92.8 
-92.9 



PBE BLYP 
Functional 




PBE BLYP 
Functional 



FIG. 1: BSE for various basis sets (left axis, vertical bars) and total FD energies (right axis, horizontal bars connected to 
guide the eye) vs. functional used in HF and PBE, BLYP DFT calculations for selected diatomic molecules containing first-row 
atoms. BSE stands for the basis set error in total energy estimated by difference i^Totai [BS] — i^Totai [FD] , FD stands for the 
finite-difference calculations. 
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One-electron energies corresponding to the six molecular orbitals of the CO molecule obtained from the local and 
nonlocal hybrids are presented in Table IVIIII The values for the la and 2a inner orbitals coincide for local and 
nonlocal methods within two significant figures (corresponding differences are about 2-4 % of orbital energy value). 
For the Itt and 5cr outer orbitals the local functionals AAm2-PBE0 and AAm2-BlLYP yield values quite close to 
those of the nonlocal hybrid methods (corresponding differences are about 2% for the Itt and smaller than 1% for 
the HOMO energy value). Whereas the SCa based local "hybrid" functional one-electron energies differ in second 
significant figure (corresponding differences are about 16% of orbital energy magnitude). 

The fact that the HOMO energies obtained from the AAm2-based local "hybrid" functional calculations are quite 
close to those of the original nonlocal ones is explained by correct asymptotic form (— f/r) of the AAm2 exchange 
potential, whereas SCa exchange potential has an asymptotic behavior similar to the corresponding DFT functional 
(PBE and B88 exchange potentials in our case) (see Ref»^ for details). 



TABLE V: Total energies from the finite-difference SC/AAm2-PBE0 and SC/AAm2-BlLYP local "hybrid" functional calcu- 
lations and comparison to the basis set PBEO and BILYP nonlocal hybrid functional self-consistent results (the differences 



between non-local hybrid and corresponding local "hybrid" functional results, A 
in parenthesis. All anergics are in Hartrees). 



-local — hyb 



local— "hyb' 



, are presented 



system _Esc-pbeo(A) -EAAm2-PBEo(A) 



BSE" 



-EsC-BlLYp(A) -E'AAm2-BlLYp(A) 



BSE'^ 



Ha -1.1689(-0.0002) 

BH -25.2491(-0.0004) 

FH -100.4005(-0.0009) 

HCI -460.6797(-0.0003) 

CN" -92.7807(-0.0015) 

CO -113.2363(-0.0015) 

NO+ -129.4622(-0.0015) 

N2 -109.4514(-0.0013) 

Fa -199.4165(-0.0022) 

BF -124.5866(-0.0016) 



-1.1691(0.0000) 
-25.2494(-0.0001) 
-100.4013(-0.0001) 
-460.6790(-0.0010) 
-92.7819(-0.0003) 
-113.2374(-0.0004) 
129.4632(-0.0005) 
109.4524(-0.0003) 
199.4180(-0.0007) 
124.5877(-0.0005) 



-1.1691 0.0000 -1.1703(-0.0001) 

-25.2495 0.0001 -25.2838(-0.0003) 

-100.4014 0.0002 -100.4728(-0.0010) 

-460.6800 0.0012 -460.8301(-0.0004) 

-92.7822 0.0002 -92.8608(-0.0015) 

-113.2378 0.0003 -113.3312(-0.0014) 

-129.4637 0.0003 -129.5637(-0.0014) 

-109.4527 0.0003 -109.5408(-0.0013) 

-199.4187 0.0004 -199.5644(-0.0021) 

-124.5882 0.0003 -124.6930(-0.0016) 



-1.1704( 0.0000) 
-25.2840(-0.0001) 
-100.4736(-0.0002) 
-460.8294(-0.0011) 
-92.8619(-0.0004) 
113.3322(-0.0004) 
129.5646(-0.0005) 
109.5417(-0.0004) 
199.5657(-0.0008) 
124.6940(-0.0006) 



-1.1704 
-25.2841 
-100.4738 
-460.8305 
-92.8623 
-113.3326 
-129.5651 
-109.5421 
-199.5665 
-124.6946 



0.0000 
0.0001 
0.0001 
0.0012 
0.0002 
0.0003 
0.0002 
0.0003 
0.0003 
0.0003 



MAE 


(0.0012) 


(0.0004) 


0.0003 


(0.0011) 


(0.0005) 


0.0003 



"Nonlocal hybrid functional calculations were performed using the largest basis sets presented in Figure [T] 
'Basis set error estimated by difference Epbe[BS] — Epbe[FD] is taken from Figure [T] 
■^Basis set error estimated by difference Eblyp[BS] — Eblyp[FD] is taken from Figure [T] 



TABLE VL Dipole moments from the finite-difference SC/AAm2-PBE0 and SC/AAm2-BlLYP local "hybrid" functional 
calculations and comparison to the basis set PBEO and BILYP nonlocal hybrid functional" self-consistent results. (Differences 
between non-local hybrid and corresponding local "hybrid" functional results, A = Mnon- local -hyb ~ Miocai-"hyb" , are presented 
in parenthesis). 

system SC-PBEO AAm2-PBE0 BSE" SC-BILYP AAm2-BlLYP BSE" 

BH -1.5248(-0.0555) -1.5379(-0.0424) -1.5803 -0.0002 -1.5119(-0.0548) -1.5254(-0.0413) -1.5667 -0.0002 
FH -1.7457(-0.0560) -1.8427(-f 0.0410) -1.8017 4-0.0001 -1.7503(-0.0572) -1.8485(4-0.0410) -1.8075 -0.0001 
HCI -1.0686(-0.0430) -1.1376(4-0.0260) -1.1116 -0.0012 -1.0463(-0.0438) -1.1170(+0.0269) -1.0901 -0.0011 
CN" 0.7183(-0.0888) 0.6159(4-0.0136) 0.6295 -0.0039 0.6836(-0.0923) 0.5805(-f 0.0108) 0.5913 -0.0049 
CO 0.2250(-0.1225) 0.1090(-0.0065) 0.1025 4-0.0001 0.1875(-0.1251) 0.0715(-0.0091) 0.0624 +0.0002 
N0+ -0.3135(-0.0767) -0.3673(-0.0229) -0.3902 -fO.OOOO -0.3277(-0.0782) -0.3816(-0.0243) -0.4059 -HO.OOOl 
BF 1.0719(-0.0713) 0.8967(+0.1039) 1.0006 +0.0002 1.0350(-0.0758) 0.8565(-f 0.1027) 0.9592 -f0.0003 



MAE 


(0.0734) 


(0.0366) 


0.0008 


(0.0753) 


(0.0366) 


0.0010 



"Nonlocal hybrid functional calculations of the electric moments wore performed using the largest basis sets presented in Figure [T] 
"Basis set error is estimated by values taken from Tables Ull IIVI 



IV. CONCLUSIONS 



Local "hybrid" functionals obtained by replacement of non-local Hartree-Fock exchange in conventional hybrids by 
local EEAX demonstrate the performance close to the original non-local hybrids. Expression for the total energy is 
the same in both cases (local and non-local hybrids), while the corresponding KS potential for local hybrids is fully 
local and multiplicative operator that brings back the local hybrids in the framework of density functional theory. 
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TABLE VII: Quadmpole moments from the finite-difference SC/AAm2-PBE0 and SC/AAm2-BlLYP local "hybrid" functional 
calculations and comparison to the basis set PBEO and BILYP nonlocal hybrid functional" self-consistent results. 



system 



SC-PBEO 



AAm2-PBE0 



BSE" 



SC-BILYP AAm2-BlLYP 



BSE" 



H2 0.5787 (-1-0.0282) 

BH -3.2961 (-0.0358) 

FH 2.2566 (-h0.0217) 

HCl 3.5797 (-h0.0603) 

CN" -5.0803 (-hO.2254) 

CO -2.0136 (-H0.0162) 

N0+ 0.5801 (-1-0.0421) 

N2 -1.5368 (+0.0925) 

F2 0.8994 (-0.0572) 

BE -3.3624 (-0.2031) 



0.6069 ( 0.0000) 
-3.2956 (-0.0363) 

2.3152 (-0.0369) 
3.6362 (-h0.0038) 
-5.0129 (-h0.1580) 
-2.0219 (-H0.0245) 
0.5742 (-1-0.0480) 
-1.5369 (+0.0926) 

0.9360 (-0.0938) 
-3.3712 (-0.1943) 



0.6069 +0.0001 
-3.3319 0.0000 
2.2783 -0.0002 
3.6400 +0.0016 
-4.8549 -0.0138 
-1.9974 -0.0004 
0.6222 -0.0003 
-1.4443 +0.0005 
0.8422 0.0000 
-3.5655 -0.0003 



0.5658 (+0.0291 
-3.3414 (-0.0315 
2.2507 (+0.0230 
3.5453 (+0.0622 
-5.1011 (+0.2362 
■2.0481 (+0.0205 
0.5489 (+0.0441 
-1.5675 (+0.0960 
0.8558 (-0.0563 
-3.4102 (-0.1886 



0.5949 (+0.0000) 0.5949 +0.0000 

-3.3369 (-0.0360) -3.3729 -0.0010 

2.3115 (-0.0378) 2.2737 -0.0003 

3.6031 (+0.0044) 3.6075 +0.0020 

-5.0321 (+0.1672) -4.8649 -0.0150 

-2.0538 (+0.0262) -2.0276 -0.0006 

0.5439 (+0.0491) 0.5930 -0.0002 

-1.5668 (+0.0953) -1.4715 +0.0003 

0.8964 (-0.0969) 0.7995 +0.0001 

-3.4048 (-0.1940) -3.5988 -0.0004 



MAE 



(0.0783) 



(0.0688) 



0.0017 



(0.0787 



(0.0707) 



0.0020 



"Nonlocal hybrid functional calculations of the electric moments were performed using the largest basis sets presented in Figure [T] 
''Basis set error is estimated by values taken from Tables [Till IIVI 



TABLE VIII: Orbital energies energies {-£i) of the CO molecule from the finite-difference SC/AAm2-PBE0 and SC/AAm2- 
BILYP local "hybrid" functional calculations and comparison to the basis set PBEO and BILYP nonlocal hybrid functional 
self-consistent results (all energies are in Hartrees). 

system SC-PBEO AAm2-PBE0 PBEO SC-BILYP AAm2-BlLYP BILYP 



la 


18.8566 


18.9016 


19.3144 


18.8765 


18.9216 


19.3350 


2a 


9.9989 


10.1093 


10.3442 


10.0215 


10.1318 


10.3675 


3a 


1.0767 


1.1536 


1.1976 


1.0692 


1.1469 


1.1907 


4a 


0.5191 


0.5910 


0.5996 


0.5137 


0.5862 


0.5947 


iTT 


0.4358 


0.5048 


0.4961 


0.4301 


0.4997 


0.4907 


5a 


0.3321 


0.3948 


0.3950 


0.3304 


0.3935 


0.3939 



In the present work the local "hybrid" functional results, total energies, orbital energies, dipole and quadrupole 
moments were compared to conventional hybrid functional values for selected diatomic molecules. The total energies 
were found to coincide for both type of functionals within the basis set error values. Multipole moments provide 
detailed information about spatial distribution of electron density. Due to sensitivity of these moments to the details 
of DFT calculation (functional employed, basis set, etc.), the discrepancies between the local and non-local hybrid 
functional results are much larger than corresponding basis set errors. The typical percentage error is usually less 
than 10% with a few exceptions. Hence the local "hybrid" functionals represent a computationally less expensive 
alternative to the non-local hybrids in framework of DFT and provide a local multiplicative exchange-correlation 
potential. Properties related to total and orbital energies, and density distributions obtained by local functionals 
were found in good agreement with those from the original non-local hybrids. 
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